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SECTION  I 
INTRODUCTION 

The  concept  of  electromagnetic  implosion  of  cylindrical  plasma  shells  (i.e., 
imploding  liners  or  hollow  Z-pinches)  to  obtain  pulsed,  high-temperature,  high- 
density  plasmas  has  been  discussed  by  Turchi  and  Baker  (ref.  1).  This  concept 
has  been  investigated  experimentally  and  theoretically  for  the  past  several  years 
in  the  SHIVA  X-ray  source  development  program  at  the  Air  Force  Weapons  Laboratory. 

The  SHIVA  imploding  liner  concept  is  distinct  from  other  imploding  liner  efforts 
(refs.  2,  3,  4)  in  that  the  idea  is  to  produce  a hot,  dense  plasma  from  the  self- 
collapsed  plasma  liner  itself,  rather  than  to  compress  a magnetic  field  preformed 
plasma  mixture. 

The  motivation  for  developing  the  electromagnetic  implosion  technique  is  to  I i 

produce  an  intense  pulsed  X-ray  source  for  laboratory  studies  of  X-ray  matter 
interaction.  The  X-ray  emission  powers  of  ~1013  watts  are  believed  possible 
with  existing  pulsed  power  facilities,  and  powers  ~1012  watts  have  already 
been  achieved  (ref.  5),  as  will  be  discussed  in  this  paper. 

A schematic  representation  of  the  basic  SHIVA  electromagnetic  implosion 
concept  is  shown  in  figure  1.  The  basic  geometry  is  that  of  a plasma  Z-pinch. 

The  load  consists  of  a large  radius,  thin  foil  cylinder  which  is  placed  in  a 
vacuum.  A fast,  high-energy  capacitor  bank  is  used  to  drive  a very  high  current 
axially  through  the  cylindrical  load.  Early  in  the  current  pulse,  the  joule 
heating  vaporizes  the  foil  and  transforms  it  into  a few  ev  temperature  plasma. 

The  JxB  force  which  results  from  the  interaction  of  the  current  with  its  self- 
magnetic  field  accelerates  the  foil-plasma  to  a high  implosion  velocity.  During 
this  acceleration  phase,  electrical  energy  from  the  capacitor  bank  is  converted 
into  kinetic  energy  of  the  imploding  plasma.  By  proper  matching  of  the  load 
and  electrical  source,  this  process  can  occur  with  good  efficiency  (20  to  40 
percent).  When  the  plasma  sheath  collides  with  itself  at  the  center  of  the 
implosion,  its  kinetic  energy  is  rapidly  converted  to  thermal  energy,  creating 
a hot,  dense  plasma  which  radiates  much  of  its  energy  in  a short  burst  of  X rays. 

To  create  a high-temperature,  high-density  plasma,  large  amounts  of  energy 
must  be  efficiently  coupled  to  the  plasma.  It  is  considerably  easier  to  deliver 
energy  to  kinetic  energy  of  a moving  plasma  than  it  is  to  directly  joule  heat  a 


Electromagnetic  Implosion  Concept.  A capacitor  bank  discharge  through  a cylindrical  foil  liner 
ionizes  and  implodes  the  liner.  Upon  collapse  at  the  center,  the  kinetic  energy  of  implosion  is 
rapidly  converted  into  plasma  thermal  energy. 
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plasma.  This  is  because  the  resistance  of  the  plasma  drops  as  the  temperature 
rises,  making  it  increasingly  more  difficult  to  couple  in  energy.  A moving 
plasma  shell,  however,  has  a substantial  dissipative  impedance  regardless  of 
temperature  because  its  inductance  is  increasing  with  time.  For  an  imploding 
foil-plasma  discharge,  carrying  axial  current  J,  the  rate  of  change  of  inductance 


where  h is  the  height  of  the  cylindrical  discharge,  uQ  is  the  magnetic  perme- 
ability (4  n X 10-7  H/m),r(t)  its  radius,  and  u(t)  the  implosion  speed.  This 

• 

L,  which  is  proportional  to  velocity  u(t),  acts  like  a resistance  with  the 
energy  dissipated  by  its  going  into  kinetic  energy.  Since  L is  much  greater 
than  the  plasma  resistance  most  of  the  energy  goes  into  kinetic  energy  and  very 
little  into  internal  energy  during  the  implosion  phase.  Thus,  with  a low  sheath 
temperature,  only  a small  amount  of  energy  is  lost  as  radiation  during  the 
implosion  phase. 

Since  the  implosion  load  draws  energy  from  the  electrical  source  and  stores 
it  as  kinetic  energy  during  the  implosion  phase  and  then  converts  it  to  thermal 
energy  at  the  center  of  the  implosion,  it  is  possible  to  obtain  significant 
increase  in  the  heating  power  over  the  power  delivered  by  the  electrical  source. 
For  an  electrical  source  delivering  energy  in  a one-microsecond  pulse  and  a 
thermal ization  phase  tens  of  nanoseconds  long,  the  plasma  can  be  heated  at  a 
power  tens  of  times  higher  than  the  output  of  the  electrical  energy  source. 

Thus,  high-temperature,  high-density  plasmas  can  be  created  with  state-of-the- 
art  power  supply  technology. 

Figure  2 shows  typical  load  current  and  the  radius  of  the  sheath  as  a func- 
tion of  time  for  a load  driven  by  a megajoule  capacitor  bank.  Note  that  the 
foil-plasma,  because  of  its  inertia,  stays  near  its  original  position  for  a 
considerable  portion  of  the  capacitor  bank  discharge  time.  During  this  time  it 
presents  a very  low  impedance;  hence,  the  current  can  rise  to  a high  value. 

The  plasma  sheath  is  seen  to  expand  initially  until  the  rapidly  rising  J X B 
driving  force  becomes  strong  enough  to  overcome  the  internal  plasma  pressure. 

The  plasma  sheath  is  then  compressed  and  accelerated  up  to  a high-implosion 
velocity. 


3 


Figure  2.  Current  and  Radius  Versus  Time  for  Imploding  Liner.  The  Two 
Curves  on  the  Radius  Versus  Time  Plot  are  for  the  Inner  and 
Outer  Boundaries  of  the  Plasma  Liner. 
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To  create  a hot,  dense  plasma  via  the  SHIVA  approach,  the  plasma  must  be 
accelerated  to  a very  high  velocity  so  that  a large  amount  of  (kinetic)  energy 
is  delivered  to  the  pinch,  and  that  energy  must  be  converted  to  thermal  energy 
in  a very  short  time.  Velocity  and  sheatr,  thickness  are  the  most  important 
parameters  because  the  thermal ization  time  is  approximately  equal  to  the  thick- 
ness of  the  sheath  divided  by  the  velocity.  For  the  highest  temperature,  the 
thermal ization  time  should  be  minimized,  i.e.,  the  plasma  should  be  heated  at 
the  highest  possible  power.  This  calls  for  minimum  sheath  thickness  and  maxi- 
mum velocity  consistent  with  good  energy  coupling  efficiency.  For  a given 
electrical  drivina  source,  there  is  an  optimum  velocity  range  for  high-energy 
transfer  efficiency.  If  high-implosion  velocities  are  achieved  by  severely 
cutting  the  sheath  mass,  very  little  kinetic  energy  is  delivered  and  though  the 
temperature  may  be  high,  the  total  radiation  energy  is  low.  The  required  higher 
velocities  may  be  better  obtained  by  using  higher-powered,  higher-energy  elec- 
trical sources.  One-dimensional  magnetohydrodynamic  (MHD)  calculations  predict 
a sheath  thickness  of  a few  millimeters.  Experimental  data  show  that  the  sheath 
is  thicker,  presumably  due  to  plasma  instabilities.  Control  of  these  instabil- 
ities, which  would  make  the  sheath  more  one-dimensional,  is  necessary  for  optimum 
load  performance. 

The  use  of  a freestanding  foil  is  critical  to  forming  a well-behaved  sheath 
that  can  run  long  distances  during  the  implosion.  This  technique  allows  the 
discharge  initiation  to  occur  away  from  the  insulator,  unlike  the  classical 
Z-pinch  in  which  initiation  occurs  across  the  insulator.  For  high-powered 
discharges  the  ultraviolet  (UV)  radiation  is  so  intense  it  can  turn  the 
insulator  surfaces  into  a conductor,  effectively  shorting  out  the  implosion 
load.  The  use  of  a freestanding  foil  in  the  SHIVA  geometry  allows  the  insulator 
to  be  hidden  out  of  sight  of  the  UV  radiation.  This  is  especially  very  impor- 
tant for  developing  a technique  which  can  be  scaled  up  to  very  high  energies 
and  powers. 
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"ECTION  II 

EXPERIMENT  DESCRIPTION 

The  imploding  liner  studies  presently  underway  involve  discharging  a 1.1 
MJ  capacitor  bank  through  low-inductance  imploding  liner  loads.  The  capacitor 
bank,  illustrated  in  figure  3,  is  described  in  detail  in  reference  6.  It 
consists  of  20  capacitor  bank  modules  spaced  around  a cross-shaped,  low-induc- 
tance, parallel  plate  transmission  line.  At  full  energy  the  top  half  of  each 
module  is  charged  to  + 50  kV  and  the  bottom  to  - 50  kV.  Hence,  when  the  bank 
is  discharged  a net  100  kV  is  applied  to  the  transmission  line.  This  technique 
allowed  the  desired  voltage  to  be  obtained  using  available  1.85  yf  capacitors. 
Special  low-inductance,  high-pressure  rail  gap  switches  were  developed  for  the 
bank  to  minimize  maintenance  without  compromising  performance.  The  bank  has  a 
capacitance  of  222  yf,  inductance  of  3 nH  and  stores  1.1  MJ  at  full  voltage. 

The  bank  has  been  used  to  deliver  peak  currents  of  up  to  20  MA  with  a current 

rise  time  of  slightly  over  one  ys.  The  load  is  located  in  the  center  of  the 
bank. 

The  two  liner  geometries  used  were  10.5  cm  radius,  1 cm  height  and  7 cm 
radius,  2 cm  height.  Both  aluminum  and  aluminized  plastic  cylindrical  foil 
liners  have  been  used,  with  areal  mass  density  of  135  to  200  yg/cm2  for 

aluminum  and  40  to  310  yg/cm2  for  aluminized  plastic.  Typical  discharge 

parameters  are  400  to  700  kJ  of  energy  (stored  at  60  to  80  kV),  peak  current 
~ 7 to  12  MA  current  rise  time  ~1  to  1.5  ys  and  implosion  velocity  ~15  to  20 
cm/us. 

The  diagnostics  included  Rogowski  coils,  voltage  probes  (ref.  7), 
single-loop  magnetic  probes  (ref.  8)}  Kerr  cell  shuttered  and  image  converter 
fast  photography,  and  radiation  detectors.  The  radiation  detectors  included 
filtered  metal  photocathodes  (X-ray  diodes  or  XRDs)  (ref.  9),  silicon  semi- 
conductor (PIN)  detectors  (ref.  10),  X-ray  pinhole  cameras,  and  convex  curved 
crystal  de  Broglie  spectrometers. 
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SECTION  III 
LOAD  CHAMBER  DESIGN 

The  evolution  of  the  load  chamber  design  is  illustrated  in  figure  4.  The 
earliest  design  --  Type  A --  was  a low-inductance  chamber  with  a 1 cm  gap  and 
a coaxial  section  to  protect  the  insulator  from  the  discharge  luminosity.  When 
evidence  was  obtained  that  electrode  ablation  was  a problem,  the  gap  was  increased 
to  2 cm  and  insulator  restrike  became  a problem.  The  Type  B chamber--without 
baffles--provided  better  insulator  protection  and  worked  reliably  at  the  60 
kV,  400  kJ  level.  At  the  80  kV,  700  kJ  level  baffles  were  needed  to  reduce 
the  discharge  radiation  incident  on  the  insulator  to  an  acceptable  level.  This 
reduction  is  6 to  10  orders  of  magnitude  in  the  visible  wavelengths.  Once  there 
was  evidence  that  baffles  worked,  the  coaxial  feature  of  the  chamber  design 
was  replaced  by  the  highly-baffled , lower-inductance  Type  C chamber  design. 

This  chamber  design  has  been  used  for  all  recent  experiments. 

Laser  shadowgram  studies  were  used  to  develop  the  baffle  design  and 
investigate  insulator  and  electrode  ablation  problems.  The  shadowgram  studies 
were  done  in  a rectangular  geometry  simulating  the  cylindrical  liner  geometry, 
i.e.,  in  a plasma  rail  gun  (figure  5).  The  studies  revealed  that  only  metal 
baffles  which  are  an  integral  part  of  the  electrodes  would  work.  The  ionized 
metal  baffle  ablatants  apparently  conducted  sufficient  current  to  be  magnetically 
confined  from  crossing  the  electrode  gap.  Notice  the  decreased  ablation  in 
regions  of  the  gap  shadowed  by  the  baffles  from  the  luminous  foil.  These 
studies  were  done  using  a single  fast  capacitor  bank  module  similar  to  those 
used  in  the  1.1  MJ  capacitor  bank. 

The  method  used  to  mount  and  insert  the  foil  into  the  electrode  gap  is 
illustrated  in  figure  6.  First,  the  portions  of  the  electrodes  inside  the 
initial  foil  radius  were  fastened  rigidly  to  a spacing  structure  which  included 
a 2 cm  diameter  post  through  holes  in  the  electrode  centers.  The  foil  was 
next  wrapped  around  the  rigidly  spaced  electrode  disks.  The  assembly  was 
lowered  into  the  chamber  and  fastened  to  the  rest  of  the  electrode  structure 
(as  illustrated).  The  original  spacing  structure  was  removed,  leaving  only  the 
foil  bridging  the  electrode  gap. 
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insertion 

removed. 
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SECTION  IV 

CURRENT,  VOLTAGE  DATA  AND  ANALYSIS 

The  current  diagnostics  normally  consisted  of  three  concentric  Rogowski 
coils--one  outside  the  vacuum--insulator  interface,  one  in  the  baffle  region 
(figure  4),  and  one  between  the  foil  and  the  innermost  baffle.  All  the  Rogowski 
loops  were  recessed  in  grooves  on  the  lower,  nominally  grounded  electrode  (which 
is  the  anode).  Typically,  the  calculated  and  measured  calibrations  agreed  to 
5 to  10  percent.  Usually  all  the  discharge  current  was  inside  the  innermost 
Rogowski  loop  until  about  100  ns  before  pinch  time. 

The  voltage  probe,  used  as  a capacitive  divider  or  as  a V orobe,  was 
located  just  outside  the  vacuum-insulator  interface.  It  was  calibrated  with  a 
low-voltage  signal  generator.  f 

Typical  current  and  voltage  data  exhibited  the  current  dips  and  voltage 
spikes  characteristic  of  pinches,  as  shown  in  figure  7.  The  measured  voltage 
spike  was  less  pronounced  than  that  for  other  pinch  devices  (such  as  dense 
plasma  focus)  due  to  the  location  of  the  voltage  probe.  All  but  3.2  nH  of  the 
system  inductance,  which  has  increased  considerably  at  pinch  time,  was  inside 
the  radius  at  which  the  voltage  was  measured.  Even  taking  this  into  account, 
the  observed  voltage  spike  was  not  as  large  as  expected  in  an  optimum  pinch. 

This  was  due  primarily  to  restrike  at  the  innermost  baffle,  about  100  ns  before 
pinch.  The  time  dependence  of  the  inductance,  L(t),  and  the  mean  current  radius 
r(t),  inferred  from  the  current  and  voltage  data  for  a fairly  good  implosion  are 
compared  to  one-dimensional  MHD  predictions  (using  the  1D-MHD  code  MAGPIE,  ref  11) 
in  figure  8.  The  discharge  current  and  voltage  are  consistent  with  one-dimensional 
prediction  until  the  last  100  to  200  ns  (or  the  last  1 to  2 cm)  of  the  implosion. 

The  inductance  and  mean  current  radius,  r,  were  inferred  from  the  current, 

I,  and  voltage,  V,  in  the  standard  fashion,  using  the  following  equations 


L = t f (V-IR)dt  (2) 


Inductance  and  Radius  Versus  Time.  Solid  Lines  are  Obtained  from  Current  and  Voltage  Data 
Dotted  lines  are  1D-MHD  predictions.  The  Two  Dotted  Line  Plots  for  Radius  Versus  Time  are 
for  the  Inner  and  Outer  Boundaries  of  the  Plasma  Liner. 
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SECTION  V 

FAST  PHOTOGRAPHY  DATA 

Fast  visible  light  photographs  were  taken  to  observe  the  position  and 
azimuthal  structure  of  the  luminous  sheath  as  a function  of  time.  Representa- 
tive 50  ns  exposure  Kerr  cell  photographs  which  were  taken  through  screen 
electrodes  during  the  implosion  are  shown  in  figure  9.  For  some  shots  the  main 
luminous  sheath  was  fairly  well  defined  with  a thickness  of  1 to  2 cm.  A 
comparison  of  visible  light  photographs  and  laser  shadowgraphs  indicated  that 
the  visible  sheath  is  opaque  and  that  the  luminosity  was  probably  from  regions 
close  to  the  ablating  electrode.  This  leads  to  an  ambiquitv  as  to  whether  the 
observed  light  is  from  the  imploding  plasma  shell  or  from  hot  portions  of  the 
electrode. 

Nevertheless,  when  photos  indicated  well-defined,  azimuthally  symmetric 
sheaths,  the  current,  voltage,  and  X-ray  data  also  indicated  good  implosions. 

In  these  cases  the  position  of  the  luminous  front  versus  time  agreed  well  with 
the  mean  current  radius  versus  time  and  with  one-dimensional  MHD  predictions 
as  shown  previously  in  figure  8.  The  final  implosion  velocities  inferred  from 
framing  and  streak  photographs  (examples  of  the  latter  are  shown  in  figure  10) 
ranged  from  15  to  20  cm/usec.  However,  the  correlation  between  visible  light 
photographs  and  load  performance  was  not  perfect,  since  ill-defined,  asymmetric 
luminous  sheaths  were  sometimes  observed  in  shots  for  which  current,  voltage, 
and  X-ray  data  indicated  good  load  performance.  Perhaps  on  such  shots  hot- 
electrode  ablatants  obscured  the  actual  imploding  plasma  sheath.  As  is  evident 
in  figures  9 and  10,  the  nain  sheath  was  preceded  by  a light  in  the  center 
(which  began  early  in  the  discharge)  for  foil  masses  greater  than  100  pg/cm2. 
The  current  associated  with  the  light  in  the  center,  i.e.,  the  current  inside 
the  main  luminous  sheath,  is  typically  10  to  50  KA  (<1%  of  the  peak  discharge 
current).  Occasionally,  luminous  streamers  or  spokes  between  the  main  sheath 
and  the  center  light  are  evident.  These  spokes  may  be  near  the  electrode 
surface  and  not  part  of  the  imploding  sheath. 

While  the  cause  of  the  early  center  light  is  not  definitely  known,  it  is 
believed  to  be  due  to  a pre-pinch  involving  the  ionized  residual  gas  in  the 
electrode  gap  (density  ~10" 5 Torr  prior  to  discharge)  and/or  a layer 
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Figure  9.  50  Nanosecond  Exposure  Kerr  Cell  Photos  of  Implosions 


FNTER  I IGHT 


VISIBLE  LIGHT  STREAK  PHOTOGRAPHS 

Streak  Photos  of  Implosions.  Image  is  diameter  of  perforated 
electrode,  swept  (in  image  converter  camera)  for  time  resolut 
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of  material  from  the  inner  surface  of  the  cylindrical  foil*  The  early  center 
light  is  not  evident  for  lower  mass  foils  (40  ug /cm2)  or  for  injected  gas 
cylinder  implosions  (ref.  12  ),  which  indicates  that  the  initiation  phase  plays 
an  important  role  in  the  formation  of  the  center  light.  In  experiments  to  date, 
the  center  light  appears  to  have  had  little  effect  on  the  implosion  of  the  main 
plasma  sheath. 


* Longmire,  Conrad  (private  communication) 
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SECTION  VI 
MAGNETIC  PROBE  DATA 

Magnetic  probe  data  was  taken  for  two-load  geometries  and  compared  with  two- 
dimensional  MHD  calculations  using  code  BRAHMA  (refs  13,  14).  The  first  data 
were  for  MJ  implosions  of  10.5  cm  radius,  1 cm  tall,  135  pg/cm2  aluminum  foils. 

The  calculations  assumed  solid  electrodes  while  the  experiments  employed  nearly 
solid  electrodes  with  3mm  diameter  access  holes  for  the  probes.  Magnetic  probes 
were  inserted  into  the  electrode  gap  at  radii  of  3,  5 and  7 cm,  with  axial 
positions  ranging  from  0.1  to  0.5  cm  above  the  lower  electrode.  Probes  at 
different  radii  were  azimuthal ly  offset  to  reduce  sheath-probe  interaction 
complications.  The  probes  were  single  turn,  ~1  mm2  and  protected  by  ceramic 
tubes. 

The  transient  magnetic  field  measurements  obtained  from  the  integrated  probe 
signals  are  shown  in  figure  11.  Comparison  of  the  absolute  values  of  reversed 
probes  indicated  that  noise  problems  were  not  severe  until  several  hundred 
nanoseconds  after  the  arrival  of  the  main  current  sheath.  The  later  arrival 
and  smaller  peak  field  seen  experimentally  is  consistent  with  a loss  of~30  per- 
cent of  the  discharge  current  at  the  insulator.  The  larger  precursor  or  foot  (pre- 
ceding the  main  rise  in  the  field)  seen  experimentally,  may  have  been  due  partly 
to  the  center  1 ight--prepinch  phenomena  and  partly  to  a diffuse  leading  edge 
to  the  current  sheath.  The  significant  agreement  between  experiment  and 
calculation  is  in  the  shape  of  the  field  versus  time  traces  following  the 
first  peak.  According  to  the  two-dimensional  MHD  calculations  (refs  13,  14), 
this  shape  is  due  to  the  plasma-electrode  interaction.  This  interaction  involves 
both  electrode  ablation  and  heat  transfer  from  the  imploding  sheath  to  the 
electrodes.  The  heat  transfer  causes  a pressure  nonuniformity  in  the  sheath, 
which  in  turn  causes  an  axial  mass  density  nonuniformity  and  results  in  an  m = 0 
instability  near  each  electrode. 

For  an  ideal  discharge  the  (integrated)  magnetic  probe  signals  should  be 
proportional  to  the  total  discharge  current  after  the  current  sheath  passes 
the  probe.  If  electrode  ablation  is  the  only  complication,  the  integrated 
probe  signal  should  exhibit  a maximum  with  the  subsequent  reduction  due  to 
shunting  of  the  discharge  current  by  the  ionized  ablatants.  The  more  compli- 
cated signal  seen  experimentally  and  calculationally  can  be  explained  in  terms 
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Figure  11.  Magnetic  Probe  Data  For  10  CM  Radius,  1 CM  Gap  Geometry. 
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of  axial -radial  structure  of  the  current  sheath  caused  by  the  heat  transfer 
perturbation  just  mentioned.  Idealized  signal  traces  for  the  corresponding 
current  sheath  structures  are  illustrated  in  figure  12. 

Magnetic  probe  data  was  also  taken  for  the  Type  C geometry  illustrated  in 
figure  4,  with  700  kJ  discharge  energy,  to  measure  the  radial  distribution  of  the 
current  density  as  time,  and  the  current  sheath  thickness.  Five  individual 
magnetic  probes  sensitive  to  EL  were  spaced  approximately  one  centimeter  apart 
radially.  The  probes  extended  5mm  into  the  electrode  gap  (from  the  grounded- 
anode  side).  The  current  density,  calculated  assuming  azimuthal  symmetry,  is 
shown  in  figure  13  as  a function  of  time  and  radius.  Note  that  for  the  early 
part  of  the  implosion,  the  current  sheath  thickness  is  < the  probe  spacing  (1  cm) 
and  that  as  the  implosion  continues,  a reversed  current  appears.  The  reversed 
current  implies  a secondary  current  path  at  about  4 cm  radius  (presumably  due 
to  ionized  ablatants  crossing  the  electrode  gap).  Anyway,  the  leading  edge  of 
the  current  shell  continues  to  implode  in  a well-defined  manner  carrying  the 
bulk  of  the  current. 

If  a secondary  current  path  at  the  4 cm  radius  has  occurred,  the  mean  current 
radius  obtained  from  current  and  voltage  data  remains  somewhat  larger  than  the 
radius  of  the  leading  edge  of  the  current  sheath  as  indicated  by  the  magnetic 
probe  data. 
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MAGNETIC  PROBE  DATA 


Current  Distribution  Versus  Time.  Obtained  from  magnetic 
probe  data  for  7 cm  radius,  2 cm  gap  geometry. 


SECTION  VII 


2D-MHD  CALCULATIONS 

Two  dimensional  calculations  of  the  implosion  were  done  using  the  single 
fluid,  axial-radial,  Eulerian  magnetohydrodynamic  computer  code  BRAHMA  (ref.  14). 
Ablation,  an  important  effect  over  time  scales  considered  here,  was  treated 
using  a simple  energy  balance  model*.  Essentially  this  model  assumes 
that  after  the  electrode  surface  reaches  its  vaporization  temperature,  all 
further  incident  energy  results  in  vaporization  of  electrode  material.  The 
reasonable  agreement  between  experimental  magnetic  probe  data  and  calculations 
suggest  that  the  treatment  of  ablation  (as  well  as  fluid  dynamics  and  magnetic 
field  transport)  is  qualitatively  accurate. 

In  both  experiments  and  calculations,  for  the  10.5  cm  radius,  1 cm  gap 
geometry,  the  electrode  ablatants  caused  considerable  (50-90%)  shunting  of  the 
discharge  current  before  the  sheath  reaches  the  center.  Recent  work  used  2 cm 
electrode  gaps,  for  which  current  shunting  through  ablatants  was  calculated  to 
be  less  than  10%  (refs.  13,  14). 

BRAHMA  calculations  predicted  significantly  better  performance  for  7 cm 
radius,  2 cm  tall  cylindrical  liner  implosions  than  for  the  10.5  cm  radius, 

1 cm  tall  liners.  The  ablation  and  heat  transfer  problems  still  occur,  but 
with  much  less  severe  consequences.  This  is  illustrated  in  figure  14,  which 
depicts  the  liner  kinetic  energy  versus  time  as  calculated  using  the  1D-MHD 
code  MAGPIE  and  the  2D-MHD  code  BRAHMA.  The  important  parameters  are  the  peak 
kinetic  energy  and  the  thermal ization  time.  Two-dimensional  effects  degrade 
performance  with  respect  to  1D-MHD  predictions,  but  power  multiplication  still 
occurs  for  the  7 cm  radius  liner. 

For  7 cm  radius,  2 cm  gap  liner  implosions,  BRAHMA  predicts  peak  densities 
~1022  cm-3,  peak  temperatures  ~150  eV,  and  an  external  effective  temperature 
~40  eV. 


★ 


Faehl , R.J.,  (private  communication) 


degree  of  power  multiplication. 
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SECTION  VIII 
X-RAY  DATA 

Low-energy  X-ray  data  was  taken  in  order  to  investigate  the  thermal ized 
plasma.  Only  axial  diagnostic  access  was  available  due  to  the  highly  baffled 
electrode  geometry. 

The  spectral  response  functions  of  bare  X-ray  photodiodes  (XRDs)  and  silicon 
semiconductor  detectors  (PINs)  are  shown  in  figure  15.  Typical  filtered  XRD 
responses  are  shown  in  figure  16.  Arrays  of  detectors  of  this  type  give  only 
approximate  spectral  information  unless  some  features  of  the  spectral  shape 
are  already  known.  Nominal  square  response  windows  used  for  such  approximate 
analysis  are  listed  in  table  1 for  the  more  frequently  used  XRD.  The  signal 
ratios  versus  blackbody  temperature  (for  assumed  blackbody  spectral  shape)  for 
these  XRDs  were  calculated  using  code  XRADAT*  and  are  shown  in  figure  17. 

Table  2 gives  the  results  for  nominal  and  blackbody  assumption  interpreta- 
tion of  such  an  array  of  XRD  signals  for  an  implosion  of  an  aluminized  (140)3) 

Kimfoil  (C16H1403)  liner  with  areal  mass  density  310  ug/cm2.  The  signal  ratios 
for  the  lower  energy  detectors  are  consistent  with  a blackbody  temperature  of 
about  38  eV. 

The  detector  signals  were  generally  consistent  with  a low-energy  continuum 
qualitatively  similar  to  a blackbody,  with  a continuum  maximum  at  ~80  to  150  eV 
photon  energy.  On  some  shots,  the  detector  signals  were  ,n  approximate  agreement 
with  a single  blackbody  temperature. 

Examples  of  the  time  history  of  some  lower  and  higher  energy  XRD  signals 
taken  for  700  kJ  discharges  are  shown  in  figure  18.  The  FWHM  is  generally 
< 100  ns.  The  X-ray  emission  powers  obtained  using  the  lower  energy  XRDs  are 
typically  ~0.5  to  1.0  x 1012  W,  assuming  isotropic  emission.  Since  the  electrical 
delivery  power  was~0.5  x 1012  W,  this  constituted  a radiation  to  electrical 
power  ration  (or  power  multiplication)  of  1 to  2.  The  yields  ranged  from  50  to 
100  kJ(total)  and  the  blackbody  temperature  ranged  from  30  to  50  eV.  In  general, 
the  highest  yielding  shots  for  700  kJ  discharge  energy  were  those  using  7 cm 
radius,  aluminum  coated  300  ug/cm2  Kimfoil  liners. 


* Johnson,  D.J.,  (private  communication) 
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Figure  16.  Filtered  X-Ray  Diode  Response  Functions.  Cathode  and 
filter  materials  identified  in  figure.  For  example, 
Ctl/FV  means  Cu  cathode  and  Formvar  anode.  FV  is  75 
ug/cm2  C5H7O2.  KF  (kimfoil)  is  250  ug/cm2  CisHiuO*. 


X-RAY  DIODE  RESPONSE  WINDOWS 
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X-RAY  DIODE  DATA  FOR  REPRESENTATIVE  SHOT 
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While  the  low-energy  X radiation  was  approximately  consistent  with  a 30  to 
50  eV  blackbody,  the  higher  energy  radiation  (photon  energy  > 1 keV)  consisted 
largely  of  A1  line  and  recombination  radiation.  The  A111+(1S2  - ISnp)  and  Al12+ 
(IS-np)  series  were  consistently  observed  using  Potassium  Acid  Phthalate  (KAP) 
and  Rubidium  Acid  Phthalate  (RAP)  bent  crystal  spectrometers.  A densitometer 
trace  with  line  identifications  is  shown  in  figure  19.  The  recombination 
continuum  is  characterized  by  temperatures  as  high  as  400  eV  (ref.  18).  The 
X-ray  yields  observed  above  1 kev  ranged  from  100  to  200  J for  aluminized 
plastic  liners  to  1 to  4 kJ  for  200  pg/cm2  aluminum  liners. 

Typical  PIN  signals,  which  are  generally  saturated,  are  shown  in  figure  20. 

The  PINs,  used  in  a helium  atmosphere  and  filtered  by  a 4 mil  Be  vacuum  window, 
were  sensitive  to  > 2 keV  X rays. 

Some  lower  energy  bent  crystal  data  (obtained  using  Lead  Stearate  bent 
crystal  with  1/8-inch  radius  of  curvature)  is  shown  in  figure  21.  The  higher 
order  lines  are  deleted.  The  Al11+  (IS2  - lS2p)  resonance  and  intercombination 
lines  have  been  seen  in  orders  2,  3,  4,  5,  6 and  7 (in  an  A1  liner  implosion). 

The  07+  (1S-2P)  line  has  been  seen  in  orders  1,  2,  3 and  4 (in  plastic  liner 
implosions).  Two  lines  (believed  to  be  first  order)  at  energies  234  and  254 
eV  have  been  consistently  observed  in  both  A1  and  aluminized  plastic  liner 
implosions.  Tentative  identifications  for  these  lines  are  Al10+  (3d-2p)  and 
A1 1 0+(3p-2s) . Except  for  these  two  lines,  the  radiation  spectrum  from  150  to 
280  eV  appears  to  be  dominated  by  continuum.  Comparison  of  the  absolute  yield 
for  photon  energies  from  150  to  280  eV  obtained  from  Lead  Stearate  data  with 
that  obtained  using  an  aluminum  cathode,  2p  thick  Kimfoil  (Cx 6Hi  403  ) filtered 
XRD  gives  factor  of  two  agreement.  The  disagreement  may  be  due  in  part  to  dif- 
ferences in  individual  crystal  reflectivities  (factor  of  2).  The  degree  of 
agreement  gives  some  confidence  in  the  use  of  XRDs  at  low  photon  energy. 

Some  representative  X-ray  pinhole  camera  photos  are  shown  in  figure  22. 

The  pinholes  were  50  to  100  p in  diameter  (measured  by  Fraunhoffer  diffraction). 
The  lowest  energy  cameras  used  to  date  used  two  layers  of  aluminized  (140  A) 

2 p thick  Kimfoil  (ClgHi403)  filters  and  Kodak  RAR2490  film,  with  sensitivity 
windows  ~200  to  280  eV  and  _>  500  eV.  Sometimes  the  radiating  region  (as  viewed 
axially)  was  donut  shaped,  with  an  outer  diameter  somewhat  less  than  the  electrode 
hole  diameter  (2  cm)  and  an  inner  diameter  of  less  than  1 cm.  Sometimes  the 
radiating  region  was  donut  shaped  with  a separate  central  hot  spot,  and 
sometimes  there  was  only  a central  hot  spot.  The  smallest  radiating  sources 
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observed  to  date  have  had  a diameter  5mm.  The  differences  in  spatial  distribu- 
tion are  not  obviously  correlated  with  other  X-ray  measurements  (yield,  power,  or 
temperature) . 

The  comparison  of  experimental  observations  and  theoretical  (2D-MHD)  predic- 
tions* of  the  X-ray  output  for  a 7cm  radius  2cm  tall,  200  ug/cm2  liner  implosion 
are  summarized  in  table  3.  Since  the  code  prediction  of  radiation  yield  is  based 
on  the  calculated  plasma  temperature  distribution  it  cannot  predict  the  experi- 
mentally observed  Al11+  and  Al12+  line  radiation.  The  temperature  contours  for 
the  "hottest"  time  in  the  pinch  are  shown  in  figure  23.  The  occurrence  of  A1 
line  radiation  is  consistent  with  the  > 100  eV  temperature  predicted  for  the 
center  of  the  pinch.  The  40  eV  contour  is  at  one  optical  depth  into  the  pinch, 
and  the  predicted  temperature  as  seen  from  outside  the  pinch  is  in  approximate 
agreement  with  low-energy  XRD  signals. 


* Roderick,  N.F.  and  T.W.  Hussey  (private  communication) 
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2D-MHD  CODE  PREDICTIONS  VERSUS  EXPERIMENTAL  RESULTS 
FOR  X-RAY  YIELD  AL  LINER  IMPLOSION 
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2D  MHO  CODE  PLASMA  TEMPERATURE  CONTOURS 


Figure  23.  Radiation  Temperature  Contours  as  Predicted  by  2D-MHD.  The  BRAHMA  predictions 
of  radiation  temperature  contours  at  the  hottest  time  in  the  pinch  are  shown 
for  one-half  of  the  electrode  gap.  The  predictions  are  for  a 200  ug/cm2  A1 
liner  implosion.  The  hot  core  temperature  predicted  by  BRAHMA  is  consistent 
with  the  occurrence  of  the  experimentally  observed  Al1*+  and  Al12+  lines. 


— 
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SECTION  IX 
DISCUSSION 

The  experimental  results  indicate  that  fairly  high-energy  liner  implosions 
have  been  achieved,  resulting  in  50  to  100  kJ,  100  ns  X-ray  bursts,  which  is 
consistent  with  2D-MHD  (BRAHMA)  predictions  Optical  and  magnetic  probe  data 
agree  qualitatively  with  BRAHMA  predictions  on  the  radial -axial  structure  of 
the  imploding  liner  sheath.  The  agreement  between  BRAHMA  predictions  and 
experimental  results  suggest  that  BRAHMA  is  a very  useful  design  tool. 

The  observed  radiation  was  consistent  with  an  30  to  50  eV  quasi  blackbody 
with  additional  higher  temperature  and  optically  thin  components.  The  experimental 
radiation  powers  of  0.5  to  1.0  x 1012W  assuming  isotropic  emission,  constitute 
a power  multiplication  of  1 to  2 over  electrical  input. 

The  power  multiplication  can,  in  principle,  be  improved  by  achieving 
sharper,  more  symmetric  imploding  liner  sheaths  in  the  final  stage  of  the 
implosion  and/or  by  using  larger  radius,  lower  mass  liners.  The  degree  of 
improvement  conceivable  is  indicated  in  table  4.  The  table  entries  were 
calculated  using  a circuit  solver  code  (DELORES)*,  with  the  assumed 
sheath  thickness  and  current  cutoff  radii  indicated.  For  the  7 cm  radius 
liner,  the  calculations,  with  2 cm  sheath  thickness  and  current  cutoff  radius 
agree  reasonably  well  with  BRAHMA  predictions  and  experiment. 

The  possible  problems  in  trying  larger  radius,  lower  mass  liner  experiments 
are  Rayleigh-Taylor  instabilities  (due  to  higher  acceleration)  and  practical 
difficulties  in  handling  liners  with  low  areal  mass  density.  The  instability 
problems  may  be  alleviated  by  one  of  several  sheath  stabilization  procedures 
under  consideration.  These  procedures  include  electrode  shaping  (using  an 
electrode  gap  which  varies  with  radius),  injection  of  axial  magnetic  field 
prior  to  the  implosion  discharge,  using  two  or  more  concentric  liners  (i.e, 
staging),  and  prevaporization  of  the  foil  liner  to  reduce  the  radial  density 
gradient  of  the  liner.  An  approach  which  might  alleviate  both  low  areal  mass 
density  and  sheath  stabilization  problems  is  to  use  an  injected  gas  liner,  per- 
haps with  preionization  prior  to  the  implosion  discharge. 


* Clark,  J.G.  (Private  Communication) 
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The  main  thrust  of  the  present  imploding  liner  development  effort  is  to 
achieve  sheath  stabilization  and  large  radius,  low  mass  implosions,  with  the 
hope  of  achieving  increased  performance  levels.  In  the  future  it  may  be  possible 
to  use  explosive  magnetic  flux  compression  generators  (ref.  23)  to  drive  foil 
cylinder  liners  (ref.  1).  Further  experiments  are  intended  to  investigate  the 
imploding  liner  concept  and  to  provide  scaling  information  for  the  design  of 
explosive  generator  driven  liners,  for  which  delivered  energy  — 1 08 J and  plasma 
formation  powers  ~101 5W,  are  conceivable  (ref.  1). 


I i 


AFWL-TR-77-252 


A1 

B 

Be 

BRAHMA 

C 

cm 

DELORES 

eV 

FWHM 

H 

h 

I 

J 

JxB 

KA 

KJ 

KV 

KAP 

MA 

MAGPIE 

MHD 

MJ 

ns 

0 

PIN 


ABBREVIATIONS,  ACRONYMS  AND  SYMBOLS 
Aluminum 

Magnetic  field  intensity 
Beryllium 

A two  dimensional  MHD  Computer  Code 

Carbon 

Centimeter 

An  imploding  liner  circuit  computer  code 
Electron  volt 

Full  width  at  half  maximum 
Hydrogen 

Electrode  gap  (height) 

Discharge  current 
Current  density 

Vector  product  of  current  density  and  magnetic  field 

Kilo-ampere 

Kilo-joule 

Kilo-volt 

Potassium  acid  Phthalate 
Mega-ampere 

A one-dimensional  MHD  computer  code 

Magneto  hydro  dynamic 

Mega-joule 

Nanosecond 

Oxygen 

Positive-Intrinsic-Negative  semi-conductor  junction 
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Electrical  resistance  of  discharge 
Rubidium  acid  phthalate 

Imploding  liner  radius  as  function  of  time  t 

Name  of  imploding  liner  program  at  Air  Force  Weapons  Laboratory 

time 

Implosion  velocity  as  function  of  time  t 

Ultraviolet 

Discharge  voltage 

Computer  code  which  calculator  detector  response  functions 

X-Ray  diode 

Angstrom 

dl/dt  (time  rate  of  change  of  current  I) 
dl/dt  (time  rate  of  change  of  inductance  L) 

Micron 

Magnetic  permeability 

Microsecond 

pi 
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